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Charles E. 

ABSTRACT 

The dynamic nuclear polarization of hydrogen nuclei ly the solid, 
effect in single crystals of samarium d.oped lanthanum magnesium nitrate 
(Sm:IMN) has been studied theoretically and experimentally. The equa- 
tions of evolution governing the dynamic nuclear polarization by the 
solid effect have been derived in detail using the spin temperature 
theory and the complete expression for the steady- state enhancement of 
the nuclear polarization has been calculated. For well- resolved solid 
effect transitions at microwave frequencies u) » cUg + the expression 
for the steady-state enhancement differs from the expression obtained by 
the rate equation approach by small terms which become zero at 
u) = + 03 ^. Experimental enhancements of the proton polarization 

were obtained for eight crystals at 9*2 (2Iz and liquid helium tempera- 
tures. The samarium concentration ranged from 0.1 percent to 1.1 per- 
cent as determined by X-ray fluorescence. A peak enhancement of l8l 
was measured for a 1.1 percent Sm:LMN crystal at 5*0° K. The maximum 
enhancements extrapolated with the theory using the experimental data 
for peak enhancement versus microwave power and correcting for leakage, 
agree with the ideal enhancement (240 in this experiment) within 
experimental error for three of the crystals. Die calculated satellite 



separation was within 6 percent of the measured separation for each of 
the enhancement curves and the peak positive and negative enhancements 
were equal for all hut two of the crystals. The nuclear spin- lattice 
relaxation time was meas\ired for one of the crystals between 1.6 K and 
K. To account for nuclear spin- lattice relaxation, spin diffusion 
theory in the rapid diffusion limit was incorporated into the results 
of the spin tenperature theory of the solid effect. The experimental 
results indicate that the spin temperature theory is a quantitatively 
correct approach for the description of dynamic nuclear polarization 
hy the solid effect for well- resolved solid effect transitions. 
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I. INTRODUCTION 


In 1953^ Overhauser^^^ proposed that the polarization of nuclei in 

a metallic solid could be changed from its thermal equilibrium value 

through the contact hyperfine interaction of the nuclei with the 

conduction electrons by saturating the electron spin resonance transition. 

The idea of dynamic polarization of nuclei was soon extended to non- 

metallic materials. Experimental verification of dynamic nuclear 

T 

polarization in metals was reported in 1956. 

Dynamic nuclear polarization (DNP) in ionic solids can occur if 
paramagnetic nuclei interact with nearby paramagnetic atoms through the 
magnetic dipole-dipole interaction. This is called DNP by the "solid 
effect" and is the subject of this thesis. Two approaches have been 
used to describe the dynamics of the "solid effect : the rate equation 
approach of Jeffries and co-workers®"^^ and of Borghini,^^ and the spin 
ten^ierature theory. ^5-17 The rate equation approach assumes that the 
"solid effect" may be described completely by following the time 
evolution of the populations of the electronic and nuclear Zeeman levels. 
This approach can lead to erroneous results if the spins are strongly 
Interacting. The spin temperature theory includes the effect of the 
strong spin-spin coupling as well as the Zeeman interactions and the 
lattice in a relatively simple way. Witliin certain limits, to be 
discussed later, the Zeeman and the spin-spin interactions and the lattice 
can each be assigned a thermodynamic reservoir having a temperature and 
specific heat. The reservoirs exchange energy with each other through 


1 
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the applied radio frequency fields, certain dipolar interactions and 

spin-lattice interactions. The expression for the steady state nuclear 

polarization resulting from the spin temperature approach reduces to the 

expression obtained by the rate equation of the Zeeman populations 

(except for a small correction term) in the limit where the electronic 

resonance linewidth is smaller than the nuclear resonance frequency. 

A material exhibiting this limit is diamagnetic lanthanum magnesium 

nitrate in which a small percentage of the lanthanum ions are replaced 

by paramagnetic samarium ions. At about 9 kOe, 9*2 GHz, and liquid 

heliijm temperatures, the electron spin resonance linewidth of samarium 

l8 

in this crystal is of the order of 5 MHz and the nuclear resonance 
frequency of the protons in the waters of hydration is 58 MHz. This 
thesis will compare the results of measurements of the dynamic polari- 
zation of protons in a number of single crystals of samarium ions in 
lanthanum magnesium nitrate (SmrLMN) with the results of the spin 
temperature theory of the "solid effect." 



II. THEORY 


A. Basic Assumptions 

A derivation of the equations governing the DNP by the "solid effect 

will be made In this section. The density matrix and spjn temperature 

theory will be the approach used In this derivation as opposed to the 

rate equation approach. The spin temperature theory has the 

advantage that the effects of the dipole -dipole Interactions may be 

Included, and yields analytical expressions for the llneshapes. A 

detailed derivation of the "solid effect" using spin temperature theory 

Is felt necessary as (l) It Is not available In the literature, (2) the 

published eqrmtlons governing the "solid effect" are either Incomplete 

or contain sign errors, and (5) all oY "ths assumptions made are not 

clearly stated. (Note In particular ref. 15-17*) 

Consider a solid containing N paramagnetic Ions per unit volume 

(which will henceforth be called "electrons ) having an effective spin 

S = 1/2 and n nuclei per unit voliome having a spin I = l/2. The 

effective Hamiltonian for this paramagnetic system when placed In a 

— > 

homogeneous magnetic field which has a steady component Hq and a high 
frequency component ^2_(^) given as^^ 

^ ■'■^IZ ■‘■^SS ■'■^II ■'■^IS ^SL ^IL ■'■^rf 

where the terms are the electron Zeeman, the nuclear Zeeman, the 
electron dipole -dipole, the nuclear dipole -dipole, the electron-nuclear 
dlpole-dlpole, the electron spin-lattice, the nuclear spin-lattice and 
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the radio frequency interaction Hamiltonians. Explicitly, they are 


5fsz 


N 

j=l 




( 2 ) 


^IZ = >'n^ 

^=l 


( 3 ) 


KlI ■ E Vl’ [4 • I; - ■ %l>] ‘5) 


k>Z=l : 


N n 


5flS " 


1 I-? [?, ■ f. 

k=l 1=1 


" ^kz 




( 6 ) 


N 


n 


^rf = ^e* Y ' ?Z 


k=l 


L. 

Z=1 


(7) 


where y (7 ) is the nuclear (electron) gyroraagnetic ratio, both 
n ® 

assumed to be inherently negative; fi is Planck' s constant divided by 

2rt' ~T is the vector distance between the S^-th and S.-th electron; 
^ ij , 

R, , is the vector distance between the Ik"^^ ^Z"^^ nuclear spin, 

K 4 

is the vector distance between the electron and the Iz“'^k 
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nuclear spin; and is the high frequency magnetic field intei*acting 

with the electrons and nuclei. The and terns are discussed 

8 

by Jeffries. 

The following assumptions are now introduced. 

1. The magnetic field is much larger than the magnetic field 

experienced by each dipole due to the neighboring dipoles. The latter 
field is called the local magnetic field. 

2. The magnitude of the electron gyromagnetic ratio 7 ^ is much 
larger than that of the nuclear gyromagnetic ratio, y^. 

5. The electron and nuclear spin-lattice relaxation times, T^ 
and Tj^, are very long conpaned to the electron and nuclear spin-spin 
relaxation times, and T^„, respectively, and there exists times 

t such that 


’■2e'’’2n « * << (| '^T ’ 

where the terms W"^ will be discussed later. 

4. The effects of the lattice may be ignored for times t. 
satisfying the ineqixality (8) and the rate of energy exchange between 
the various spin systems and the hi^ frequency magnetic field H 2 (t) 
calculated. The effect of the lattice is then reintroduced. 

5. The electron and nuclear resonance lines are homogeneously 

20 

broadened. 

Consider the large constant magnetic field to be applied along 


the z-axis in the laboratory reference frame. Then 


6 


X 


sz 


7 

e 


H 


N 

Z_J zj 
j=l 


(9) 


and 



= 7 -ft 
n o 



( 10 ) 


Define 


o> 

e 


7 H 
e o 



= 


Z=1 


( 11 ) 

( 12 ) 

(15) 


(iM 


where cUg and are the electron and nuclear Larmor frequencies, 

and S and I are the z-conroonents of the total spin of S- and 
z z 

I-spins, respectively. Since 7 g 7n> thus if an 

oscillating field %(t) of frequency o) is of the order tOg, then the 

direct interaction of the nuclear spins with this field is negligible. 

— > 

We apply this oscillating field in , a plane perpendicular to H^, so 
that equation (7) becomes 


7 


N 

^rf = * % Z (^xk ^ Syk 

k=l 

or 

3 C _ = 7 -h H-, (S cos cut + S sin cut) , 

ri G X X j 

where and are defined in a fashion similar to in equation 

(15). Defining ciij^ = and noting the identity 

g-icifizt g g+iofizt ^ g cot + sin cut , 

X ^ y 



where 

AkZ = ^zk Szi (1 - 5 cos^ \i) 

®kz = ■ I 5 ^kzHsJs'^ + S^S'J) 

^kZ = ■ 2 ®kZ ®kZ ® (^zk^Z + ^zZ^k^ 

®kZ = • I ®kZ ®kZ (^zk^Z ^zz\) 
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where 0^^^ 
the z-axis 

JC and 3f can he similarly written. The effective Hamiltonian now 
IS II 

has the form 

-iofi t iccSjjt 

5f ' = + ^SS ^IS ^II ■*' 

where the prime indicates that the spin-lattice interactions have been 
omitted. Note that the Hamiltonian has been written in frequency units, 
i.e. , fi = 1. 

The time evolution of the system described by the Hamiltonian can 

be determined from the equation of motion of the density matrix, p, 

. 21 

for the system. The density matrix obeys the relation 

to . i[p, W] , ( 18 ) 

again taking h = 1. 

Following Redfield,^^ we transform the density matrix by a unitary 
transformation, R(t), to a frame of reference rotating about the 


E. 




5 sln^ e,, ^..3 




z 2 ” “ 

^kz = - r ®ki ® 


(l6b) 


and are the polar coordinates of the vector 


is parallel to H^^ and 


S, = S , ± i S, 
k xk ky 
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laboratory z-axis at a frequency o). Thus 

p^(t) = R(t) p(t) R"^(t) , 


where P (t) 


is the density matrix in the rotating frame and 


R(t) 


^ioBzt 


( 20 ) 


We do not transform the nuclear spins because they are not subject 
to saturating rf fields in these experiments. See Wollan ^ for a 
discussion of this point. 

The eqiiation of motion for the density matrix then becomes 

( 21 ) 


where 


?C (t) 




(D I 
n z 


{■ 


+ COj^S^ +R(t) 3fgg +5fjg 


} 


R'^(t) 


( 22 ) 


and 


A = iUg - as . 


( 25 ) 




is not altered by the unitary transformation since 



JfgS and 3Cis are of the form 
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R(t) 3fgs R“^(t) 

= JCgg + 

}]^ss 

m 


R(t) 5Cjg R-^(t) 


m 

(24) 

where the s\ims are over the values 

m = ±1, 

SS 

and X° are the 
IS 


secular parts of the dipolar Hamiltonians, i.e., those which commute 
with R(t): 




= 0 , 


(25) 


and dff™ and are the nonsecalar parts of the dipolar interactions. 

SS IS 

The nonsecular terms in equation oscillate at frequencies ±2o) 

in the rotating reference frame and may he ignored as discussed hy 
Redfield^^ and Goldman. Then from equations (l6) and (25) 



2 

y 

e 


N 


k>Z=l V 


\ 



(26) 


and 


N n 


X 


IS 


- ^ y 

/_, L, \ k 


kZ^zk^zI ^kZ^zk^Z ^ ^kZ®zk^ 


k=l Z=1 


0 


(27) 


where 




(28) 
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= - t ^e>'n «kl 


0,,, cos 0 J 5 .J e 


■l'>lti ^ -5 


( 29 ) 


and the star indicates complex conjugate. The Hamiltonian in the 
rdtating frame is reduced to 


3C^ = ASg + %Iz + ^OlSx +3fgS + ^IS 


o _ p 


( 50 ) 


We define 




( 51 ) 


and 


V = ^II * 


( 52 ) 


Then 


3f " = 3C^ + V . 


( 55 ) 


Transforming to the interaction representation hy a vmitary trans- 
formation 

,R/ 


p^(t) = u-Ht) p^(t) U(t) 


( 54 ) 


where 


U(t) me ^ , 


( 55 ) 


one has 


Ir - {■='’ 


( 56 ) 
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where 

V^(t) = u"^(t)V U (t) . (57) 


Eqviation (36) may be integrated and then iterated to obtain the form 



and this resiilt differentiated to obtain 


= ip(o), V^(tj] . ^ ... jT 

n=l 

X p.. [p"(o), vi(t„)], vi(Vi)]. • < 


I, » R 

Note that p (o) = P • 

The basis of the spin temperatxire theory is to assume at this point 
that the density matrix in the rotating frame will have the canonical 
form for times t satisfying equation (8). The justification for this 
assumption is discussed in references 16, 22, and 2h. For this system, 
it is assumed that the commuting terms AS^, and 5Cgg can be 

considered as thermodynamic reservoirs having inverse temperatures a, 

3, and 7, respectively. Each reservoir is considered to reach internal 
equilibrium in a time of the order of its spin-spin relaxation time. 
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Therefore, we assume that 


-a(t)ASz - P{t)o:hlz " ?'(t)3fsS 

P^(t) = 


The inverse temperature, a, is defined as 


-a(t)AS„ - P(t)a:^l2 


( 1 + 0 ) 



where k is the Boltzmann constant and Tg^ is the electron Zeeman 
spin temperature. The inverse temperatures P and y are similarly- 
defined. 

The electron-nuclear dipole-dipole reservoir and the nuclear 
dipole-dipole reservoir have been ignored in equation (l+O) as -rfell as 
the off-diagonal part of the density matrix. The electron-nuclear 
dipole-dipole reservoir (i.e., the equation (27)) 

contributes to the inhomogeneous broadening of the resonance lines and 
is assumed to be negligible. This term will also give rise to a 
barrier to spin diffusion and this will be treated phenomenologically 
when spin diffusion is taken into account in section II F. The off- 
diagonal terms in are retained in V as a perturbation. 

The nuclear dipole-dipole interaction may be ignored since its 
frequency (-'^0 kHz) is very small compared to those of the electron 
Zeeman, nuclear Zeeman, and electron dipole -dipole interactions in the 
rotating frame (all -5 MHz). This term will be important only when spin 


Ik 


diffusion is taken into account. The size of the contribution made by 

x6 

the off-diagonal part of the density matrix is discussed by Goldman 
and for the time scale considered here, the off-diagonal part may be 
ignored. With these approximations, we note that 

P^(t) = P^(t) . 

The representation used for the electron system is one in which 
5^g and are diagonal. Since the nuclear dipole-dipole interaction 

has been neglected, the nuclei are considered independent particles and 
the representation used is one in which all the single particle operators 
1-^ and are diagonal. 

In the high temperature approximation, one has 

p^(t) = T]"^ |l - a(t)AS2, - 3(t) - 7(t)dfCggJ , (4l) 

where t) is the normalizing factor. This approximation is necessary to 
make the calc\jlation tractable, and corresponds to the experimental 
conditions encountered in this work. 

Note that the electron and nuclear magnetizations along the z-axis 
are proportional to the ensemble average of the operators and I^, 

respectively. Using the density matrix, these values are found as 

<S2> = Trp^(t)S2 

<I > = Trp’^(t)! 

Z z. 


and 
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where < > indicates the canonical ensemble average value and Tr 
indicates a sum over the diagonal elements in a space that includes both 
the electron and nuclear eigenstates. It follows that the time rate of 
chajige of these expectation values are (in the high ten^jerature limit) 


d<S > 
z 

dt 


= Tr 


dp^(t) 

dt 


Sz = 


' ^ dt 


( 1 + 2 ) 


and 


d^ 

dt 


= Tr 


dP^(t) 

“dT“ 


-%T^Trl/ 


dt 


i^3) 


It can be noted from (4-2) and (45) that the changes in the magnetization 
are proportional to the changes of the corresponding inverse temperatures 
of the reservoirs. The calciolatio;i of the time evolution of the inverse 
temperatures using equations (39)> (4l), (42), and (45) follows. 


l6 


5 , Derivation of the Equation of Motion for the 
Electron Zeeman Inverse Spin 
Temperature , a( t ) 

From equations (39) and. (^2), 


-(ATrs/) ^ 
2 dt 


i TiTr[pP, V(t)]]s2 + (i)^ TTr j' dt'j^P^, V(t’)]^ V(t^ 

0 ^ 


+ (: 


0 0 
t 


t' r 

• fdt" 

J 

rp, v(t")], v(t 

0 

t' 

L 

t” 

rs 


• /dt" j 

’ df" 


0 0 


- 

)], v(t") 

1 , V(t') 

, V(t) 


- 



S + ♦ • • 

z 




where the superscript I, indicating the interaction representation, has 
been dropped. Ass\aming the trace and the time integrals commute, the 
following terms must he evaluated! 


= t,Tt [p^ V(t)]s^ 



Using the identity 
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Tt[a,^C = Tr a|b,c] , (^6) 

one has 

Cl = T]Tr [s^, P^]v(t) = 0 , (*+7) 

since 



Using the same identity 


Cg = 


= t)Tt , v(t')]j |[v(t), sjj . 

p^(t) can he written in the following way 

p^ = |l - (a - y) - O - - rjfeJ . 

Recalling that 


(1^8) 


{^9) 


v(t) = ^S^(t) + 3Cjg(t) , 


where 


S (t) and 
x^y 




are in the interaction representation, we 


have 


Cg = A(a - y) Tr Sy(t') Sy(t) 

+ i a>^.u„ (3 - 7) Tr S^(t) (t')] 

+ i /c^iTr Sy (t)[3C^, V(t')] . (50) 
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The second term is zero as can easily be seen using eq\iation (46) and 
noting that 


Cs (t), ij = 0 . 


Using 


1 Tr S^(t) [JC„, V(tO = 


and integrating noted in (44), we get 

t t 

P dVC2 = -^3^^A(a - 7) J ~ 

0 0 

t 

0 

The second term in (51) is also zero. This is shown by letting 
T = t - t ' ; then 

J dt' ^ Tr[Sy(t - t')v] = -J 

0 0 


But 


Tr Sy(r)3C^g = Tr Sy(r) 


r N n 

I I<\z * 4i 

k=l i=l 


= 0 (55) 


since it has no diagonal elements. Noting the identity 
6 



19 


Sy.(T) = - S^(T) Sin AT + Sy(T) cos AT , 




where 


13(?„ T 

Sx,y(’-) ■ - ^x,y - 


-‘4 


( 55 ) 


the remaining term in (5l) 'becomes 
t 


0 t 

= d.T ^ [rrS^Sy(T) cos At - TtS^S^(t) sin At] . 


Note that 


TrS^Sy(T) = 0 , 

since a rotation of it about the y-axis sends into and leaves 

Sy(T) unchanged. Since the trace is invariant the result follows. Then 

* 


The term TrS^S^(T) is proportional to the correlation function of the 
transverse magnetization^ with a correlation time Since t is 

assumed to be very much larger than Ige^ the upper limit of the inte- 
gration may be taken to infinity. Thus 


pOO 

-ax^J dr ^ [TrSxSy(T)] = a^[rrSx^(T) sin At] 


T=oo 


= 0 , 


T=0 
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since the correlation function is zero at t = oo. Equation (5l) is 
reduced to 


dt* Cg = 


A(a - 7) / dt’ 


CO 

/' 


TrSySyCt 


- t') 


Using exactly the same arguments as used above 



dt' Cp = 


00 

-A(a - 7) j' dr TrSy^(x) 


cos At , 


(56) 


C may be calculated using some of the relations developed in 

3 

finding C^. Expanding we get 

Cj = -co^^A(a - 7) TrSy(t")[3f°g(t’), Sy(t)] 

+ i^^^^ (p - 7) Tr [l^, 5Cis(t")][kis(^' Sy(t)] 

+ i7^^ Tr v(t")][v(t’), Sy(t)] + 0 . (57) 


The first term in (57) can be written as 


TrSy(t")[jC°g(t'), Sy(t)] = Tr[Sy(t), Sy(t")] 3C jg(t ’ ) . 

This term is zero for the same reason that equation (53) zero. The 
second term in (57) can be written as 

" Tr^x(^’)^ Sy(t)], 3Cis(t") = 0 
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since 

j[S,(t' ), Sy(t)], = 0 . 

The third term in eqxiation (57) is easily seen to he zero hy 

iitS 

applying the unitary operator e ^ to the argument of the tirace. 

Since S^ commutes with and 3Cjg(t), the rotation sends ®y(^) 

into -S (t) with the res\ilt that 

y 

The fourth term in equation (57) when integrated over t" becomes 

irCDlJ' dt" Trji^,V(t"^|j(tO,Sy(t^ = Tr|v(t') - V(0)^ [v(t'),Sy(t)] 

= Tr[y(t’), V(t’)] Sy(t) 

- 70 i^ Tr V(t’)[Sy(t), V] (59) 

The first term above is obviously zero, and the second term can be 
expanded to give 

Tr V(t’)[Sy(t), v3 = Tr sJK^(V) + [Sj;(t’), 

+ 01^ Tr Sy(t)[»^, 3tis(t')] + 0 • 

The first term on the ri^t hand side is zero for the same reason that 
eQ^uation (53) zero and the second term is zero using the identical 
argument leading to eqxaation (58)* Therefore, to order of = 0. 
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C can be evaluated using arguments identical with those used in 

4 

evaluating C, , C^, and C , along with the condition that the terms of 

1 3 

the order of and and higher are neglected. The result is 

(a - 7) Tr[Sy(f"), 3fyt")][3C^(V), Sy(t)] 

+ (u . r)rTrfe,3fis(t"')], S^(t"^[<s(V), Sy(t)] 

* 0- > > ®y<‘0) <60 


It is now assumed that each nuclear spin interacts with only one 
electron. For dilute systems, this is an excellent approximation. Thus 
the sample is broken up into N equivalent "spheres of influence, 
each containing one electron and ^ protons. In this approximation, 
one has 





n 



I. 

J 


(61) 


where the summation is over one sphere of influence, and is for the 

jth pyoton referenced to an electron at the center of the sphere. 
Integrating as indicated in equation (44) yields 
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X t t' t" 

Jc^it = (a - r) Tr ^ /«" 

J=1 000 

|^i%(t"-t’) -i03„ (t" - t')J 


X S^(t"’ ) S^(t) 


n 

N 


i- i 


+ (3 


I A*' 


X Sy(t") S^(t) 


[.‘ 


J=1 000 

io^ (t'” - t') -la^ (t"’ - t’)J 


(62) 


- e 


We define I (o^) as 
t 

I(%) 


^ ^ ' icjj (t” - V ) 

" ' dt'" S^(t"’ ) Sjj(t) 


dt" e 


t" 


0 0 


and integrate this expression by parts 


i(%) 


io) J 


t ( ioi (t" - t’) rt" 

J dt'” S^(t'” ) S^(t) 


dt ' < e 


n 0 I 0 

t' 


t"=t' 

t"=0 


,/dt" 


io^ (t" - V) 


S^(t") S^(t)> , 


which reduces to 



2k 


«“n) - 

^ 0 


f 


(t” - 1 ’ ) 


/ r i^v, . 

dt' J dt" e S^(t") S^(t) . 


0 0 


Integrating the second term on the right in the same manner, 

t t' 

f dt' f dt” S (t") S (t) 
n icn J J ^ ^ 


0 0 


. p i03 (f - t) 

/ dt’ e Sj^(t’) Sx(t) 


+ P dt' S (f) Sjj,(t) . 


The first term on the right of equation (62) is 
n 

(a - r) Tr + l(-»„)) - 


j=l 


n 

N t 

m • < 

2 


2£o) ^ 2 + - r 


»2 
n 


J=1 

n 

N 


iTiT / dt' S^(t' - t)S^ cos A(t' - t) 


^ ^ I k 


"iti: 

J 0 


n 


J=1 


X / dt* S (f - t) S^fcoB (cD + A)(t' - t) + cos (tu - A)(t' - t) 

-V X x|_ n J 


0 
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We use the same procedure to evaluate the second terra on the right 
of equation (62) as was used to obtain equation (65), and assume that 


e . 

J 


may be replaced by an average over the sphere of influence of 


the electron, i.e., an e such that 


n 

N 


. 1=1 


( 61 +) 


In section II F, we evaluate for the conditions corresponding to 

our experiments. Equation (62) becomes 


J Ci^dt = -A (a - 7) I (A) (Tr s/) 

+ A (a - 7) |[w"('\ - A) + W*-(oj^ + Tr 
- 0)^ (p - 7) |[w"(“n - A) - W+(“n + A)j Tr (65) 


where 


(JU. 


W°(A) = -^ g (A) , 


(66) 


W (o)^ ± A) 


- 

= ^ 4- g (“ ± A) 




2 2 n 


(67) 


The electron spin resonance absorption line shape function g(^*^) is 


defined as 
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g(“) 

and normalized such that 



(68) 


00 



= 2rt . 


The fact that 


00 

Tr Ij I~j J dT ^(t) Sx cos orr ^ Tr ^ dr ^(t) cos cut 


Tr S„2 


Tr S„ 


was used to obtain equation (65) from equation (65). W° and W* are 
the transition probabilities per unit time associated with the allowed 
electron spin resonance transition and with the solid effect transitions, 
respectively. 

Combining eqcaation (44) with equations (56) and (65)- (68) leads to 


da 

dt 


. 1 _ n 2e y 

\ “nV 


+ 2. i w"(co - A) 


l^cu^P - Ax - (o^n - A)/J 
_ 2. 1 W+^O) + A) < CO p + Ax - (co + A)7? 

N A n ^ n n j 


(69) 


where 2. ^ is a correction factor to the allowed transition probability 

N 2 

per unit time due to the Cj. term in the perturbation expansion. Since 
it is usually much less than unity, we drop it from this point on. The 
terms a, p, and 7 in equation (69) are rigorously a(0), p(0), and 
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7(0). For times t satisfying equation (8), they are approximately 
equal to a(t), 3(t), and 7 (t), respectively, and we use this approxi 

mation henceforth. 


C. Derivation of the Equation of Motion for the 
Nuclear Zeeman Inverse Spin Temperatiire , p(t) 

The derivation of the equation of motion for p(t) is similar to 
the derivation of the equation of motion for a(t). Starting from 


equations (39) and (^3), we get 

t -j 

.^(Tr l/) ^ = in Trp, V(t)] + (i)^n Tr J V(t')J, V(t)J 

t t* — » 

+ (1)5,, Tr j df f dt" ^[p^ V(t”)], V(V)J, V(t) 

t t’ t" 

+ (i)\ TrJ'diV J it" J dV" 


j- 0 0 0 

r“ _ 


|ip^ v(t”')], v(t"j, v(t>) 

, V(t) 


- 


I + • • • 

^ (TO) 


We 


define traces analogous to the C's of equation (^5) 


C^' = H Tr V(tj] 

Cg' = n Tr V(t’)], V(t)J 



( 71 ) 
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' and C ^ ' 


are easily shown to be zero. Evaluating C 2 ' 


C^' = H Tr [p^ v(f )][v(t), ij 


= -(13 - 7)<»„ Tr [l^, 3C°g(f)][3Cyt), Ij 


where the other terms are zero. 

Using eqmtion (2?), carrying out the commutation relations and 
integrating C 2 ' as indicated in equation .( 70) , we get 


N n 

f dt' C 2 ' = (3 - r) Tr ^ Y, "ij4j 

i,k=l j-1 


SzAk (•") p 


-icu T 

n 


The product (t) is related to the correlation function for the 

z-component of the individual electron spins^ which involves a correlation 
time of the order of T 2 g- The limits on the integral may then be taken 
as infinity. This expression reduces to 


J dt' C 2 ' = (3 - 7)(Tr 


(72) 


where 


n N 


-icu T 


Tr X Z ^zAk ^ 


W 


IN j=l i,k=3^ 
2 n 


Tr S. 


(73) 
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We now use the sphere of influence model to deduce an estimate for 
equation (75). In this limit only i = k terms will contribute 
significantly to equation ('73) and one gets 

N 


\ 


Tr ^ J dT (t) e 

2 « ^ .00 

w 

1 2 ^ 2 


■ICD^T 


Tr S„ 


where £ was 


^ defined in equation (64). We now assume the approximation 
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r “ ... "^^n^ r °° ~ \ 

Tr 2^ j dT (T) e - Tr / dr (t) e 

i=l -00 


Using equation (68) we then have 


u » £_ e(o) ) 
1 2 n'^ 


(74) 


If one assumes a Lorentzian lineshape,^^ equation (74) leads to the 
same result as reference 25. W^ as defined by equation (75) is the 
transition probability per unit time resulting from the fluctuation of 
the z -component of the electron spins due to the electron dipole -dipole 
(3("_) interaction. The effect that this term has on the nuclear 

OD 

relaxation has been discussed in reference 25. 

C4' is found using the same arguments used to obtain eq\aation (60). 
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Cl,' 


= i(a - y) Tr ^|^Sy(t"' ) ) ^ 

+ [[Sy(f"), 3c5s(t")], S^(f)][3(53(t), ijj 


- CJD cu^ (3 - 7) Tr 
n 1 


[V K?s(t'"0. 


Q(Sg(t), 1^3 


Making the same assumptions which led to equation (61), and 

integrating Ci^' as indicated by eqviation (70) yields 

n 

t t' t" 

j dt C^' = iA(a - 7) ^ 2-1 I j 

J=1 




0 0 


/ i% 

X S^(t"')Sy(t*) 

n 

N 

2 m,, V 


iQ^(t"-t) -io) 
- e 




t" 


>n(P - y) ^ Tr ^ 

J=1 


2 1"!:!: /"df r dt" rdt"' 

J oj J J 

0 0 0 

/icon(t"'-t) -ia^(t"’-t)\ , ^ 

Sy(t")Sy(f )(e^ +e j. (75) 


Integrating hy parts as done before, and using equations (64), (67), 
and (68), eqiiation (75) becomes 


I 


dt 


Cj^' = -A(Tr 1^) (W" - W+) (a - 7 ) + ^^^(Tr 1^) (w" + W^) (3 "^^)^* 


Combining equations (72) and (76) with equation (70) yields 
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dt 


D. Derivation of the Ecmation of Motion for the 
Electron Dipole-Dipole Inverse 
Temperatvire , 7 ( t ) 

The derivation of the equation of motion for the electron dipole- 
dipole inverse temperatiire, 7 (t), can be accomplished using the fact 
that energy is conserved in the rotating reference frame for times t 
satisfying equation ( 8 ), provided that 

<K’> » <V> , 
o 

where 3C and V are defined in equations (3l) aril (52)^ respectively, 
o 

Then 


= -W^(3 - 7) - i W - A) I u^3 - Ax - (u)^ - A)W 

- -L w'*' (ui + A) < rn 3 + Ax - (<n + A)7? • 
n n “ J 


( 77 ) 


d<5t> 


dt 


= 0 = Tr(AS„ + 




O) I 
n z 




dt 




so that 



Substituting equations ( 69 ) and (77) into ( 78 ) and using 
(Tr Ig^/Tr S^,^) = n/N, we get 
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dt 


+ 2. -0- W-I (3 - 7) + W°(A) (a - 7) 

N ^c. 

+ 2. ^ W" (oi - A) (cn^3 - Ax - - A)7> 

N 01^2 J 

n ^ y+ /cn B + Ax - (^„ + A) 7/ 

N 0^2 " L ” J 


(T9) 


where 


2 


“i - 


Tr S 


The term is related to the linewidth of the electron spin resonance 

, ^ 14,26 

line. It is found to be 



where Mg is the second moment of the electron spin resonance line. 

E. Effects of the Lattice Reservoir 
The effect or the lattice vibrations (phonons) is introduced by 
assuming that each spin-lattice relaxation process takes place exponen- 
tially with a characteristic spin-lattice relaxation time, and that 
these processes may be added to the equations of evolution for a, 3^ 
and Thus equations (69), (??)> and (79) become 
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da 

dt 


-W°(A) (a - 7 ) + - - W(o) - A) 3 - A3- - - A)7 

N A " 




^ i W‘^(oj + A) 
N A ^ 
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3 + At - (o)^ + A)7 






( 80 a) 


dt 


= -W^^(3 - 7 ) - ^ - A) - At - - A)^ 


- ^ + A) 

n V 


oa^ 3 + Ax - (tii^ + A)7 


'1 -- 
J T 

n 


(3 - 3^) (80b) 


.2 . ^ CD_ 


^ = A^^ W°(A)(a - 7 ) + a ^ w (3 - 7) 
dt n\ 2 JN 






+ S. 

N 


(OD - a) 

" — W"(o)^ - A) { cu„3 - At 




( CU £ 


(O)^ - A) 7 


') 




+ At - + A )7 


i -k 


(7 - 3 J 


( 80 c) 


where Tg, (see section II F), and are the spin-lattice 

relaxation times of the electron-Zeeman, nuclear -Zeeman, and electron 
dipole-dipole reservoirs,®^ respectively, and 3 l is the 

inverse temperature of the lattice. Note that the electron-Zeeman 
inverse temperature does not approach the inverse lattice temperature. 
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■but the "colder" temperature — 3 t. This is a result of the transfor- 

A ^ 

mation to the rotating frame which affects the electron Zeeman 
interaction, but does not affect the nuclear Zeeman and the secular 
electron dipole -dipole interactions. 

F. Nuclear Spin Diffusion 

The nuclear dipole-dipole interaction, igiiored in the 

derivation of the kinetic equations governing the solid effect. This 
interaction becomes important if N « n, for there will be large 
regions in the solid where the electron -nuclear dipole-dipole interaction 
is very small compared to the nuclear dipole-dipole interaction. The 

nuclear magnetization in these regions will be spatially transported by 

+ - 

energy conserving transitions resulting from the terms ll' 

This process is called nuclear spin diffusion. Only a brief discussion 

of spin diffusion in the rapid diffusion limit will be given here and 

l4 

will closely follow the treatment given by Abragam and Borghini. 

27-51 

More elaborate discussions of spin diffusion are given elsewhere. 

Each electron can be considered to interact only with the nuclei 
within its sphere of influence having a radius R defined by 

E . (!^ , . (81) 

where N is the n\imber of paramagnetic ions per unit volume. An 
important parameter is the pseudopotential radius b, which is the 
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distance from the electron at which a nucleus has approximately the same 
probability of being relaxed by the electron as of being flipped by the 

28 

JCj-j terms due to neighboring nuclei, and is given by 


b = 0.68 




(82) 


The term C is 


5 2 2,2 

C = — 7 7 1.1 


2 m 2 " 


10 'e 'n 1 + ^ 2 T 
n e 


(85) 


and D is the diffusion constant given approximately by 
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1 a 
D = -i- 


10 T, 


2n 


(81) 


Here a is the average distance between nuclei and is the nuclear 

spin-spin relaxation time. Another important parameter in spin 
diffusion theory is the diffusion barrier radius bo, defined^^ as the 
distance from the paramagnetic ion at which the magnetic field at the 
site of the nuclear spin due to the ion is equal to the nuclear line- 
width and given very approximately as 




The regime of 


"rapid diffusion" 


28,50 


( 85 ) 


defined by the inequality 
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a, b « « R , (86) 

corresponds to the conditions for the crystals and temperat\ires reported 
in this thesis. 

Detailed calculations^^^^^^^^^^^ show that in the limit of rapid 
diffusion, equation (80b) can be written as 

^ = -W (p - Pq) (8T) 

dt 


where 


W 


+ W 




A) + W (ol - A) + ^ 

m 


T, 


n 


(88) 


and 


Po 




^ / j (89) 


The bars correspond to angulsir and radial averages within a spherical 
shell whose radius r is bounded by 


bo < r < R . 

2 

In this limit and within these approximations, 6 , defined in equation 
(64), becomes 
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2 2 
>'e ’'n 



> 


( 90 ) 


and is the appropriate expression to use in equations (6?)^ (7^)^ and 
all three equations (80). Henceforth, we assume that this is done, and 
omit the bars over and W“, but not over One can show 


that 


14 , 30,31 


A. = iil + — . (91) 

T 5 "b 5 
n 

The first term is due to nuclear relaxation by the desired paramagnetic 
impurity (the samarium ions, in our case), and the second term is 
"leakage" relaxation by other impurities and by other modes of spin- 
lattice interaction. The physical picture is that the relaxation and 
rf transitions caused by W tend to drive the inverse temperature 3 
to a steady state value given by 3^ same time the spin 

diffusion mechanism attempts to keep the total nuclear magnetization 
spatially homogeneous. 

It is left to determine the measured spin-lattice relaxation rate. 
In the absence of any radio frequency fields, equation ( 87 ) can 
be written as 


^ = -W, (3 - 7 ) - f =-'\ (P - Pl^ • 

['^J 


The time evolution of 


7 


must be taken into account also, i.e.. 


38 


2 

. p) . - p ) . 

dt Td 


(95) 


The general solution for 3 exhibits two time constants. There are 
several limiting cases for which the longer time constant can be easily 
evaluated; this will be the measured T^- One limit is 


d’ 


in which case 




A second limit is the situation in which 


t 


V Uvj " Vui 


... f ■ < 3-7 

In this case 7 will come to quasi-equilibrium ^i.e., 
t such that 


O) in times 


Thus, equation (93) may be solved for 7, which is then substituted 
into equation (92). The latter can be rewritten as 




dt T, 


n 



59 


where 


W. 


T. 


nOD, 


n 


1 + Tj) 


If 


1 » 


ncu. 
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^NcDl‘ 


2 )\ 


equation (9^+) becomes 


1 _ ^ 1 _ 

T T 
n 


+ . 


On the other hand, if 


one obtains 




Tn nO)^2 Tp 


Combining equations (91) (9^)^ 'w® Set 


(9M 


(95) 


( 96 ) 
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1_ 

Tn 


kn nc ^ 1 ^ 

3 "^nl 


Wi 


ruo, 


1 + 


n 




2 


( 97 ) 


which is consistent with the crystals and temperatures considered in 
this thesis. The last term on the right side of equations (9^) through 
(97) is a result of the coupling of the nuclear spins with the electron 

dipole-dipole reservoir.^5 note that l/Tj) is often taJcen^*^ to 

he 2/Tg. 



G, Theoretical D^mamic Nuclear Polarization R esults 
The dynamic nuclear polarization, usually expressed in terns of 
the steady- state enchancement E defined as 


E = 

can be deduced from equations (80) and (8 t) - (8?) by setting all the 
time derivatives equal to aero and solving for p. The cumbersome 
general expression can be deduced from Cramer's Rule and can be found 
in appendix A. We find E for several limiting cases. 

Consider first the case of well- resolved solid effect transitions, 
that is, W- , 5^ 0,- W^ , W" = 0. One finds for ± A) « 0 the 

result 


1 ^ W^T 


E 


o)n 


n 


(98) 


n Tg 


1 + W-T„ 1 + - 


N T 


n 


where is the measTired nuclear spin-lattice relaxation time given 

in eauation (97)- For saturating microwave powers » Di 

equation (98) reduces to 


E = 


- ^e 

+ — 

T 

1 + S^ 


(99) 


N T, 


n 


The tem - — has been called the leakage factor in the litera- 
ture.^^^^5 3? this term is small compared to unity, then equation (99) 


becomes 



k2 

E = ± ^ . (100) 

The nuclei have in effect been given the thermal equilibrium temperature 
and polarization of the eleotrona. mis first limit, to which equation (98) 
applies, corresponds to the experiments reported in this thesis. 

A second limit is that of unresolved solid effect transitions, 

ignoring all leakage factor effects. In this case, one can deduce 

[ ’ 

E — ~ 3“ 

(l + W3_Tn + W"Tn+wX) 

This is the Abragam and Borghinil^ result with the addition of the 
factors. At low microwave powers, this becomes 

(jDe (W “ 

E « 1 + ^ 

1 + WiTn 

whereas, for saturating levels of microwave power, one gets 



««e 
E « — 


u> A 
n 




The leakage factor corrections to these equations requires the full 
solution to equations (80) and ( 87 ) - ( 89 ) and is given in appendix A. 


III. TOfPKRIMEMAL APPARATUS AMD PR0CEDI3BE 


Figure 1 iB a picture of the experimental apparatus used in this 
study, me experimental equipment consisted of a nuclear magnetic 
resonance (WE) spectrometer, an electron spin resonance (ESE) spectre- 
meter, an electromagnet, and a cryogenic system. 

A. Nuclear Magnetic Resonance Spe ctrometer 
The NMR spectrometer is shown schematically in figure 2. The 1®® 
detector was of the "Q-meter" type. A detailed discussion of the theory 
of operation and of the limitations of this detector is given in 
reference 32. The rf oscillator vas a constant voltage device, and a 
ramp generator was used to sweep the frequency. A "line-stretcher" 
was used to maintain a half wavelength line between the sample coil and 
the variable capacitor that made up the tank circuit of the Q-meter. 

The sample coil consisted of two turns of teflon coated number 56 
copper wire wound on the crystal. The coil leads were fed through a 
small hole drilled in the end wall of the microwave cavity and connected 
to the coaxial cable. Care was taken to have the plane of the sample 
coil parallel to the microwave magnetic field to prevent coupling of 
the microwave power out of the cavity through the NMR spectrometer. 
Standard lock-in detection was used to record the NMR spectra. 

B, Electron Spin Resonance Spectrometer 
Tte ESR spectrometer used in this experiment is shown schematically 
in figure 3. The spectrometer was a standard circuit employing 



mtmni 



■r 

-p" 


Figure 1.- Apparatus used in the DNP experiments. 
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Figure 2.- Block diagram of the nuclear magnetic resonance spectrometer 
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COILS 


Figure 3-“ Block diagram of the electron spin resonance spectrometer. 
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a three-port circulator to utilize all the power furnished hy the 

klystron. The klystron was a Varian Model V 58 reflex klystron rated at 

500 milliwatts at X-hand. The klystron frequency was locked to the 

resonance frequency of the sample cavity using an automatic frequency 

55 

control circTiit similar to the Berry and Benton circuit. A 
sketch of the cryogenic probe containing the sample cavity and variable 
coupler are shown in figure 4. The sample cavity was machined from 
brass and gold plated and resonated in the cylindrical TE 111 mode. A 
wire was soldered on a diameter of the sample cavity to short one of the 
two degenerate resonant modes. A small cylindrical section of teflon was 
screwed to an end wall of the cavity and used as a crystal mount. The 
variable coupler used was similar to that discussed in reference 54. 

The saii 5 )le cavity and variable coupler were attached to a length of thin 
wall stainless steel type 504 waveguide^ and this was soldered to a 
flange which bolted to the top of the cryostat. Standard lock-in 
detection was used to record the ESR spectra. 

C. Magnet 

The magnet employed was a Varian 25 cm electromagnet having a 
6.7 cm air gap, with Fieldial^^^ and Hall probe control. Coils were 
positioned around the pole pieces to provide the magnetic field 
modulation necessary for lock-in detection of the signals. 

D. Cryogenic System 

The apparatus used to meas\ire and control the teniperat\ire of the 
liquid helium bath is shown in figure 5- Th.e temperat\are of the helium 




Figure 4.- Sketch of the cryogenic sample probe. 
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Figure 5.- Schematic of the cryogenic system. 
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bath was monitored by measuring the vapor pressure of the liquid helixim 
and converting the vapor pressure to temperature using the 195^ Helium-^)- 
Temperature Scale. Temperatures below 4.2° K were maintained in the 
8 liter heliiim dewar by pumping over the liquid with mechanical roughing 
pumps having a total p\imping capacity of 10 liters/sec. A constant 
vapor pressure was maintained by a. manostat connected between the helium 
bath and the roughing pump as shown. 

E. Sample Preparation 

The samples were prepared starting with the rare earth oxides 
La203 and SnigOj obtained from the Lindsay Division of American Potash 
and Chemical Corporation. Both the 1*205 and Sm205 had a 99-99 percent 
rare earth purity. Both oxides contained the naturally occurring 
isotopic compositions. The rare earth nitrates were synthesized separately 
using reagent grade concentrated nitric acid in the reactions j 

(Sm^I»a)203 + 6HNO3 -♦ 2(Sm,La) (NOj)^ + 5H2O . 

The pure rare earth nitrates were then added to reagent grade magnesium 
nitrate in the reactions: 

2(La,Sm) (N03)5 + 5Mg(N05)2 + 24H20 -» (La,Sm)2 Mg5(N03)^2 * 24H20 . 

The pure an and La double nitrate solutions were allowed to saturate at 
room temperature and then mixed to yield aqueous solutions of samarium 
doped lanthanum magnesium nitrate (Sm:LMN). The doped solutions were 
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poured into beakers, placed into a desiccator having concentrated 
sulfuric acid acting as a desiccant, and maintained at 0° C. 

Single crystals of the rare earth double nitrates weighing approxi- 
mately 200 mg. could be grown in about 50 hours. The growth habit is a 
flat hexagonal plate with the crystal symmetry axis perpendicixLar to the 
plate. The detailed crystal structure of these salts can be found 
elsewhere. ^ 

F. Experimental Procedure 

Small slits were cut into the edges of the sample crystals to 
align the two t\jrn NMR coil. The crystal was glued to the teflon 
sample mount on the end wall of the microwave cavity and the probe 
immersed in liquid helium. The samari\am ESR line was located and 
recorded. The angle 6 between the crystal symmetry axis and the 
dc magnetic field was determined and fovind to be approximately 7° in 
all but one case. 

With the microwave power off, a thermal equilibrium NMR signal was 
recorded by sweeping frequency at a given temperature and magnetic 
field. Thermal equilibrium nuclear polarizations at different temper- 
atures and magnetic fields were calculated using this result and the 
Brillouin function. With the dc magnetic field approximately 50 Oe 
below the main ESR line the microwave power was turned on and a dynamic 
equilibri\mi allowed to bd established between the various reservoirs. 
The enhanced NMR signal was recorded by sweeping frequency. The dc 
magnetic field was changed and dfter a few nuclear spin-lattice. 
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relaxation times the NMR spectrum again recorded. We measured the NMR 
derivative peak-to-peak height using magnetic field modulation with a 
peak-to-peak amplitude of 2.5 Oe and a frequency of 50 Hz. Because of 
signal-to-noise problems, the thermal equilibrium NMR measurements had 
a relative error of ±20 percent, while enhanced NMR signals had a 

relative error of ±5 - 5 percent. 

The nuclear spin-lattice relaxation times were measured by 
enhancing the NMR signal with microwave power, then turning the micro- 
wave power off and recording the decay of the NMR derivative peak. 

Any frequency drift in the rf oscillator during the measurements of 

was corrected by monitoring the frequency and ai>plying a correction 
voltage to the voltage- controlled oscillator. The data was plotted on 
semi-log paper and determined from the slope. 



IV. EXPERDffiNTAL RESULTS 


The experimental results for eight samarium doped lanthanum 
magnesium nitrate crystals are summarized in table I. The positions 
of the hyperfine lines of and were determinedly. 57 for the 

applied magnetic fields and frequencies listed in table I and found to 
be far removed from the solid effect and main ESR transitions, and 
thus they should have no effect on the DNP process except possibly 

through leakage relaxation Tj^^* 

The magnetic field Hq reported as item 5 in the table is the 
magnetic field value halfway between the enhancement extremums and is 
assumed to be the field at which the main ESR transition occurs for 
the frequency reported as item U. The value of Ho reported is within 
+5 Oe of the value measured when the ESR transition was located and 
recorded. The variation of the position of the main ESR line is due 

to the backlash in the field control. 

The angle, 0, between the ciystal symmetry axis and the applied 

magnetic field was determined using the values 

= 0.756 +0.005 
= 0.565 +0.10 

reported in reference 58 . The value for g„ was measured for a number 
of crystals and found to have a value of 0.752* The value of g^ 
could not be measured due to the field strength limitations of the 
magnet. The parameter g of the electronic system is defined by 
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TABLE I. Simmary of experimental parameters: Item 1 la the crystal designation; item 2 Is the measured tej5)erature In degrees Kelvin; Item 5 Is the magnetic field 

in oersteds halfway between enhancement extrema; Item !*• is the microwave frequency in gigahertz used to induce the forbidden transitions; item 5 is the angle in 
degrees between the applied magnetic field Hq, and the crystal symmetry axis; Item 6 is the electron spin-lattice relaxation time in seconds taken from reference 59} 
item 7 is the measured nuclear spin-lattice relaxation time unless otherwise noted; item 8 (9) is the measured derivative peak-to-peak line width In oersteds of the 
ESR (NMR) transition; item 10 is the satellite separation in oersteds, measured and calcuiAted; item 11 (12) is the maximum positive (negative) enhancement measured; 
item 15 gives the percent dilution of samarium in the crystal solution and the measured dilution in the crystal; item 14 is the leakage factor calculated from items 
6 7, and 15; item 15 is the infinite power enhancement, calculated from the dependence of enhancement on microwave power, neglecting the leakage factor; item 16 is 
discussed In the text; item 17 is the enhancement corrected for leakage; and item l8 is the ideal enhancement and is discussed in the text. 


1 . 

CRYSTAL 

1 

2 



4 


5 


6 



7 


8 

2 . 

TEMPEHATORE (°K) 

1.23 

3.06 

1.81 

1.87 

1.22 

1.22 

3.1 

1.85 

1.22 

3.0 

1.21 

3.19 

2.6 

1.21 

3 . 

Ho (Oe) 

8915 

8925 

8903 

8900 

8952 

8926 

8906 

8898 

8908 

8872 

8935 

8918 

8879 

8912 

L. 

Ve (GHz) 

— 

9.112 

9.J30 

9.121 

9.155 

9.156 

9.128 

9.118 

9.120 

9.089 

9.117 

9.127 

9.081 

9.110 

5 . 

0 (degrees) 

— 

0 

7 

7 

7 

6 

6 

6 

7 

7 

7 

7 

7 

8 

6. 

Te (see) 

2.5 X 10-5 

0.0018 

0.10 

0.10 

2.5 X 10*5 

2.5 X 10-5 

0.0015 

0.10 

2.5 X 10-5 

0.002 

2.5 X 10*5 

3.1 X 10"** 

0.011 

2.5 X 10-5 

7 . 

Tn(see) 

12 

100a 

250 a 

200 

15 

15“ 

100® 

250® 

15® 

108® 

15 

57 

151 

16 

8. 

OH^p (Oe) 

3.5 

5.2 

.... 

5 

5 

5 

5 

5 

5 

5 

6 

6 

6 

1.6 

9 . 

OHgp (Oe) 

.... 

.... 

.... 

15 

10 

12.5 

10 

— 

11 

9.5 

8.1 


8.7 

.... 

10. 

SATELLITE SEPARATIOB (Oe) 
Measured 
Calculated 

73 

71 
71. 1 


71 

71.2 

71 

71.6 

72 

71 

71 

71 

71 

71 

71 

71 

75 

71 

71 

7‘*.3 

71 

71.3 

75 

71.5 

78.1 

71.2 

11. 

*^max 

61 

78 

22 

38 

36 

25 

50 

17 

56 

117 

30 

61 

95 

55 

12. 


61 

80 

.... 

35 

10 

21 

50 

17 

56 

181 

32 

79 

122 

35 

15 . 

SAMARIUM DILUTION (*«) 
Solution 

20 

15 

15 

20 

20 

25 

25 

25 

25 

25 

20 

20 

20 

20 


Measured^ 



0.1 

0.3 

0.5 

0.8 

0.8 

0.8 

1.1 

1.1 

0.3 

0.5 

0.3 


11. 

(1 + n Te/HTn) 

— 

— 

10.6 

5 

1.01 

1.01 

1.05 

2.2 

1.01 

l.Ol 

1.01 

1.01 

1.6 

— 

15 . 

E' 

— 

117 

25.8 


197 

— 

77 

— 

187 

202 

100 

115 

— 

— 

16. 

So 

— 

0.11 

0.070 


1.21 

— 

0.52 

— 

2.10 

0.122 

2.15 

0.81 

— 

— 

17 . 

E' (1 + n Tg/NTn) 

— 

— 

215 

190 

199 

.... 

81 

— 

189 

210 

101 

120 

— 

— 

18. 

^ideal 

210 

210 

2I0 

210 

210 

2lC 

2l0 

210 

2I0 

2I0 

210 

210 

210 

210 


STaken from figure 11. 

'^asured samarium concentrations were made by 
X-ray flotirescence. 
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vdiere 3^^ is the Bohr magneton. For axially symmetric crystals like 

8 19 

Sm:LMN, one can show ’ ^ 

-)l/2 

cos^6 + sin^6 

There was indication from the literature^^ that the samarium ion 

is highly rejected by the host lattice of IMN» A single crystal of 

samarium magnesium nitrate was grown and used as a standard to determine 

the actual concentration in the an:IMN using X-ray fluorescence. The 

solution and measured dilutions of an:LMN for each of the crystals is 

listed as item 13 in table I. The measured dilution of samarium in IMN 

T 

was used to calculate the leakage factor — — — ^ which occurs in the 

JN 

denominator of equation ( 98 ). The electron spin-lattice relaxation 
time Tg was calculated from the expression 

-i- = 3.1J. T + 1.5 X 10"^ T^ + 1.6 X (101) 

Tg 

which was determined experimentally by Larson and Jeffries^^ for 
dilute an:IMT in the parallel orientation under conditions similar 
to ours. The crystals we used were not all in the parallel orienta- 
tion when the data was taken, but Tg for SmrlMT does not vary 

< o 

appreciably in this temperature range for angles 0 ^ 10 . 
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A reproduction of the derivative of an ESR absorption signal is 
given in figure 6. Most of the ciystals used in this experiment 
erfiibited a similar asymmetry. The derivative peak-to-peaJt values 

AH® of the ESR absorption line for the crystals are listed in table I. 

PP 

Figures 7 and 9 are enhancement curves for two crystals. The 
enhancement is the ratio of dynamic to thermal equilibrium proton polar- 
ization. The solid lines are smooth curves through the data points. 
Figures 8 and 10 are the variations of the peak enhancements as a 
function of microwave power for the same two crystals. It is evident 
that the peak enhancement is limited by microwave power. The enhance- 
Di0nt versus relative microwave power data was fit by computer to the 
equation 

E = E’ ^ — (102) 

S+ S3 


This is equation (98) expressed in terms of E (= p/Sl) 
defined as 


E’ 


“n 


1 + — — 
N T„ 


(103) 


(The one (l) in the numerator of equation (98) is neglected compared 
to other terms in the numerator. ) S is the microwave power in 
arbitrary units (maximum experimental power corresponds to S = 1), 
and Sq is a constant determined in the computerized least squares 


fit. 
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10 Oe 




Figure 6.- The derivative of the ESR absorption line for crystal k. 
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lOOL 
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CRYSTAL 6 
T=3.0 °K 



Figure 10.- Enhancement of the polarization peaJi as a function of 
relative microwave power for crystal 6. 
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The values of E' and Sq determined by the computer fit are 
given as items 15 and l6, respectively, in table I. The curves dravn 
in figures 8 and 10 are equation (102) using the E' and values 

thus determined. 

The experimental peak enhancements can be compared to the ideal 
enhancement, which is the ratio of the electron and nuclear 

Larmour frequencies, that is. 


E 


ideal ’ 


(ICA-) 


by correcting for the leakage factor (nT^/KT^). The result of this 
correction is given as item 1? in table I. 

The nuclear spin-lattice relaxation time, T^, was measured as a 
function of bath temperature between 1.6° K and i<-.2° K at a constant 
magnetic field value of 897 O Oe using crystal 7- These data are 

shovm in figure 11. 

A theoretical enhancement curve can be drawn if a line shape 
function for the ESR absorption is assumed. Equation (98) can be 
written in the form 


E = 


1 + E’ 


ISojg(O) 


1 + 


S\ g(AH) 

SoJs(O) 


(105) 
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using the terms defined in equations (102) and (lOj). We use here 
the shorthand notation g(AH) for e{o^ t A) which appeared in 
equation (67). Equation (102) was written for the peaic enhancement 
at g(0). The curve in figure 12 is equation (105) plotted for 
crystal 4 with maximum microwave power (S = 1 in our units), assuming 
a Gaussian line- shape function 

g(AH) ^ -2(AH)2/(AH^p)2 
g(0) 

with parameters from table I for ciystal 4 . The circles are the 


measured DNP enhancements. 





V. DISCUSSION OF RESULTS AM) CONCLUSIONS 


The major objectives of this research were to derive the equations 
governing the dynamic nuclear polarization by the solid effect using 
the spin temperature theory (eqs. (8o)) and to test this theory in the 
limit of well- resolved solid effect transitions. Single ciystals of 
samarium doped lanthanum magnesium nitrate were used since the solid 
QffQcti transitions of the samarium ion in the host lattice of lanthanum 
magnesium nitrate are well resolved from the main ESR transition. 

The theory predicts a maximum, "ideal" enhancement given 

in equation (104) for complete saturation at the center of the solid 
effect transition with negligible leakage. For the situations realized 
in this thesis 

^ideal = 

Data of the peak enhancement as a function of relative microwave 
power were obtained for crystals 2 through 7* each case, the data 
indicated that insufficient power was available to completely saturate 
the solid effect transition (see, for example, figs. 8 and 10). The 
maximum enhancement can be predicted from this data if it is assianed 
that the enhancement as a function of power is given by equation (98)^ 
which was rewritten as equation (102). The extrapolated value of 
Infinite microwave power enhancement (E* in eq. (102)) is given as 
item 15 in table I. Satiiration was almost complete at 1.8l° K for 
crystal 3 since the observed peak enhancement was 22 and the value E' 
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6 ? 


was 25.8. After the leakage factor was taken into account, the extrapo- 
lated maximum enhancement (item 17 in table I) was 24-5* This value is 
well within experimental error (i'20 percent ) of the theoretical value, 
240. The extrapolated maximum enhancements of crystals 4 and 6 were 
20 percent lower than E^^^ggl corrected for leakage, which is 

the lower limit of our experimental, error. The remaining crystals were 
not analyzed as above, or showed low enhancements. The reasons for the 
latter- are not understood. 

The theory predicts that for well- resolved solid effect transitions 
the enhancement extremums shoTold be symmetrically situated about the 
Tnfl.i n ESR transition and at 

OD = o)^ i . 
e n 

Since it was difficult to establish the position of the center of the 
main ESR transition with respect to the positions of the peak enhance- 
ments, the solid effect separation (labeled as satellite separation, 
item 10 of table I) was used. This separation in oersteds should be 

This value was calculated and is shown in t^'tla I. The measured 
separation is in excellent agreement with the calculated separation 
for each crystal and is well within the 5 Oe variation due to backlash 
in the incremental field controller. . - 

The theory also predicts that for saturation conditions the abso- 
lute value of the enhancement extremums should be the same, that is. 
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Tiax ^Wix * 

A comparison of items 11 and 12 in table I indicates that this is true 
for most of the crystals and temperatures. However, in the values for 
crystal 6 at 5*0° K and crystal 7 at the two lower temperatures, 
is up to 22 percent less than E" . This behavior has been observed 
before^ and was attributed to random drifting in the microwave 
apparatus . 

The theoretical expression for the steady- state enhancement cnirve 
is a function of the ESR line- shape and microwave power. If these are 
known, a curve can be drawn predicting not only the positions of the 
enhancement extremums but also the shape of the enhancement curve. 
However, the ESR line- shape function was not computed from the measured 
ESR derivative curve in this experiment. Instead, a theoretical 
enhancement curve was drawn assuming a Gaussian line-shape for the 
ESR line and using the values of E' and Sq determined from the 
enhancement versus microwave power data for crystal k} this is shown 
in figure 12 along with the experimental enhancement data for this 
crystal. The enhancement curve determined in this way is very sensitive 
to the derivative pealt-to-peak width of the main ESR line. 

The enhancement curves drawn through the experimental data in 
figures 7 and 9 do not agree with the theory in the region between the 
enhancement extremums. At approximately 10 Oe on either side of the 
enhancement extremums, the value of the enhancement should be +1. This 
deviation is attributed to the distorted line- shape of these crystals. 
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The theory of spin diffusion was incorporated into the spin 
temperature theory of the solid effect in order to account for the 
nuclear spin-lattice relaxation time. Only the case of rapid spin 
diffusion was considered, with the restilt that the nuclear spin- lattice 
relaxation rate should be given equation (97)^ can be seen from 
the analysis in section II F and the data in table II. If we assume 
that there is no leakage relaxation ** since the first term 

on the right side of equation (97) is negligible compared to the last 
term, then for crystal 7 we find 

2 

" ■ wi ^ 4 ^’ 

where we have used the relation^** 2Ij, = T^. The nuclear spin-lattice 
relaxation time calculated from this equation results in times of the 
order of 10^ seconds at K and 10^ seconds at 1.6° K. These 

calculated relaxation times differ from the experimentally deteimined 
relaxation times by one order of magnitude at 4.2^ K and three orders 
of magnitude at 1.6° K. This deviation may be due to other para- 
magnetic impurities in the crystal that would result in a non- 
negligible leakage relaxation time. This mechanism would result in 
a measured relaxation time shorter than the theory would predict. 
InsTifficient work was done to understand this discrepancy. 

In conclusion, the equations governing the dynamic nuclear polar- 
ization by the solid effect using the spin temperature theory have 
been derived- For well- resolved solid effect transitions at microwave 
frequencies cd » cOg t the expression for the steady- state 



TABLE II.- NUCLEAR RELAXATION AND SPIN DIFFUSION PARAMETERS 

FOR Sm:LMN 


1 . Samarium 

dilution {%) 

0.1 

1-1 

2 . R (nm) 

5-2 

2.4 

5 - n/N 

2 lf ,000 

2,200 , 

h. T^ (sec)^ 

2 . 5 x 10"5 - 0.10 

2.5x10"^ - 0.10 

5 . Tgg (sec)'^ 

5.6x10"® 

3.6x10"® 

6 . T2n (sec) 

4.5x10"^ 

4.3x10"^ 

7. D (nm^/sec) 

950 

930 

8. C (nm^/sec)*^ 

7 . 2 x 10"5 - 1.8x10"^ 

7 . 2 x 10"5 - 1.8x10"' 

9. a (nm) 

0.2 

0.2 

10. b (nm) 

3.6x10"^ - 4.5 x10"5 

5.6x10"^ - 4.5x10" 

11. b (nm) 

1.2 

1.2 


12 . 

O'. 

Wq_ (sec"l) 

2 . 0 x 10 "^ 

2 . 2 x 10 "^ 

15- 

W, (sec"l) 

3.1 X 10 ^ 

5.1 X 10 ^ 





14. 

(sec"l)® 

2.9x 10"5 - 7.2x10"^ 

5.2x10"^ - 7.9x10"® 


®X-ray fluorescence. 

^rom reference 59 * 

CEstimated from ESR peak-to-peaX linewldth assuming that ESR 
line is completely . homogeneously broadened. 

%or » 1. 

®Estimated from equation ( 91 )- 


TO 
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enhancement differs from the expression obtained by the rate equation 
approach by small terms Which become zero at os = Og + o)j^. The 
maximum enhancements extrapolated with the theory from those obtained 
experimentally agree with the ideal enhancement within experimental 
error for three of the crystals. The calculated satellite separation 
was within 6 percent of the measured separation for each of the enhance' 
ment curves, and the peak positive and negative enhancements were equal 
for all but two of the crystals. These results indicate that the spin 
temperature theory is a quantitatively correct approach for the 
description of dynamic nuclear polarization by the solid effect for 
well- resolved solid effect transitions. 
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VII. APPEND K 


A. Complete Expression of the Steady- Stat e Enhancement 
The complete expression for the enhancement of the nuclear polar- 
ization can be obtained from eqviations (8o) by setting all the time 
derivatives equal to zero and solving for p using Cramer’s Rule. 

The result is 


E 


Pl Q 


■where 


>- Ae ^■^.n'^e n ~ ^ j 

-I “e “n+^,n^e n (“n ^d\ 


+ W^, 


+ w°tJi + ^ ^ (i + I T 

“l V 


“l 


? T. 


— / A CO0 Ip n Tp 

■H W°T,(WiT„ + WX 1 -T i: ^ i T f 

\ “l ® ^ 


7 ^ 
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+ W. T (l + ^ ^ ~ 
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+ W°Tg ^1 + ^ ^^(1 + W^T^ + W+T^ + W-T^) 
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